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Abstract 

Firstly, a systematic procedure is derived for obtaining three-dimensional 
bound-state equations from four-dimensional ones. Unlike "quasi-potential 
approaches" this procedure does not involve the use of delta-function con- 
straints on the relative four-momentum. In the absence of negative-energy 
states, the kernels of the three-dimensional equations derived by this technique 
may be represented as sums of time-ordered perturbation theory diagrams. 
Consequently, such equations have two major advantages over quasi-potential 
equations: they may easily be written down in any Lorentz frame, and they 
include the meson-retardation effects present in the original four-dimensional 
equation. Secondly, a simple four-dimensional equation with the correct one- 
body limit is obtained by a reorganization of the generalized ladder Bethe- 
Salpeter kernel. Thirdly, our approach to deriving three-dimensional equa- 
tions is applied to this four-dimensional equation, thus yielding a retarded 
interaction for use in the three-dimensional bound-state equation of Wallace 
and Mandelzweig. The resulting three-dimensional equation has the correct 
one-body limit and may be systematically improved upon. The quality of 
the three-dimensional equation, and our general technique for deriving such 
equations, is then tested by calculating bound-state properties in a scalar field 
theory using six different bound-state equations. It is found that equations 
obtained using the method espoused here approximate the wave functions ob- 
tained from their parent four-dimensional equations significantly better than 
the corresponding quasi-potential equations do. 

I. INTRODUCTION 

With the advent of experimental facilities such as CEBAF, which are capable of probing 
hadronic systems at energies where relativistic effects become important, the development 
of theoretical frameworks in which consistent relativistic calculations of few-body hadronic 
systems can be performed is no longer merely desirable, it is essential. 
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An obvious starting point for such a calculation in a two-body system (such as the 

deuteron) is the Bethe-Salpeter equation (BSE) for the four- dimensional, covariant, two-to- 
two amplitude T, 

T = K + KG T, (1.1) 
where Go is the free two-particle propagator, which in our convention is 

G (p' 1 ,p , 2 ;Pi,P2) = t(27r) 8 6^(p' 1 - Pl )S^( P f 2 -p 2 )d 1 ( Pl )d 2 (p 2 ), (1.2) 

with 

. / x _ A t(Pi) K(Pi) /, ~\ 

a i{Pi — o / \ , • i / N ■' \ L - 6 

where 



ei(Pi) = VP? + m i' ( L4 ) 

| „ I r , for spin-zero particles, 

* Wt) I ±£ ' (p ' ) 7 7 7 V p ' +m % for spin-half particles, 1 J 

7] is a positive infinitessimal and K is the Bethe-Salpeter kernel In principle should 

include all two-particle irreducible two-to-two Feynman graphs. Solution of (|1 . 1|) with the 
full two-particle irreducible kernel is impractical and usually resort is made to the ladder 
approximation Some undesirable features attend this approximation (see JIIJ for a 



full discussion), not the least of which is the fact that the ladder BSE does not have the 
correct one-body limit []TT|. (An equation is said to have the correct one-body limit if, when 
the mass of one particle is taken to infinity, the equation reduces to the Klein-Gordon or 
Dirac equation for the light particle moving in the static field of the now infinitely massive 
source.) By contrast, the full BSE ( |1 . 1|) does possess the correct one-body limit. 

Three-dimensional quasi-potential equations (QPEs) are realized by rewriting the Bethe- 
Salpeter equation for the two-particle amplitude ( |1.1|) as a pair of coupled equations: 

T = U + UgT } (1.6) 
U = K + K(G -g)U, (1.7) 

where the propagator g is arbitrary. If g is chosen to contain a one- dimensional delta function 
constraining the relative four-momentum, then Eq. ( |1.6| ) becomes a three-dimensional equa- 
tion for the amplitude T. Examples of this approach include the equations of Blankenbecler- 
Sugar JTJI and Gross [P] , |I3(] . In principle the driving term U should be determined by solving 



the four-dimensional integral equation ( |1.7|) . However, given the difficulty of solving such 
equations, and our lack of knowledge about the best form of K for hadronic physics, usually 
U is chosen to be a one-boson exchange interaction: 

U = Vobe- (1.8) 
The t-matrix then obeys the three-dimensional integral equation 
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T = Vqbe + Vqbe gT. 



(1.9) 



A reasonable description of few-hadron systems is obtained by fitting coupling constants 
and some mass parameters of Vqbe to the nucleon-nucleon scattering data. (Examples of 
this approach include, but are by no means limited to, Refs. [|I4| -|T6|| . ) 

The QPE formalism of Gross obeys the one-body limit 1 1 1| . Another variant of the 
quasi-potential approach which respects the one-body limit was derived by Wallace and 



Mandelzweig in References ||17H19|. We provide a generalization of that formalism in this 
paper. As originally derived, the QPE of Wallace and Mandelzweig contains, for the nucleon- 
nucleon system, a one-boson exchange potential which is instant in the center-of-mass frame. 
This involves the use of a constraint on the relative four-momentum, p, of the form 

p-P = 0, (1.10) 

where P is the total two-body four-momentum. In the center-of-mass frame, the equation 
for the deuteron vertex function takes the form: 

^WM = V[ ns tG E T^WM, (1-11) 

where Get refers to the Wallace-Mandelzweig choice for the propagator g. The formalism 
is covariant since the constraint ( |1.10|) is expressed covariantly and the pair of equations 



( |1.11| ) and (|1.7|) is equivalent to the bound-state BSE: 

T = KG T. (1.12) 



Eq. ( |1.11| ) for the cm frame deuteron vertex function has been used in recent work by 



Devine and Wallace P0| , [21|| . This formalism may be developed in a manifestly covariant 
fashion following the technique of Fuda . This involves writing all four- vectors in terms 
of their components parallel and perpendicular to the total four-momentum: 

p = P\\P + p±_ (1-13) 

where p\\ = p ■ P, with P = P j \[P^. 

A fundamental flaw exists in quasipotential formalisms that are based on a form for 
g which contains a delta function. It is generally impossible to systematically correct the 
lowest-order approximation by use of Eq. (|1.7|) because unphysical singularities arise. In 
two-body hadronic systems, such as the deuteron, calculations may still be pursued using 
equations such as (|1.11| ) which are themselves free of these singularities. On the other hand, 



in any electromagnetic reaction where the photon momentum q is non-zero, the constraint 
( 1.10| ) cannot hold for the initial and final-state relative four-momenta in the cm frame [20 



Therefore a boosted vertex function is required for electromagnetic matrix element calcu- 
lations, or, for that matter, for any QPE calculation in a three-body problem. A boost 
equation for the quasipotential may be deduced from Eq. ( |1.7|) , but unphysical singularities 
arise in the boosted interaction These are removed from the theory if the full result 



for U, as defined by (|1.7|) , is used in the quasi-potential equation, but no truncation of U at 
any finite order in G — g is free of singularities. 

Therefore, in this paper we seek a general procedure for the reduction of four- dimensional 
equations which does not involve the use of delta functions. We develop a procedure by 
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which a four-dimensional equation may be approximated by a three-dimensional equation 
and the interaction in that equation improved systematically. The technique is presented 
in Section [TT]. It has a close connection to the work of Klein [p3|-|26|| on three-dimensional 



reductions of four- dimensional equations, and to standard time-ordered perturbation theory, 
as discussed in Section fTJ. 

In general the interaction in the three-dimensional equations we discuss is not covariant. 
But, since the application of our delta-function- free reduction technique to infinite order 
produces an equation equivalent to the original four-dimensional equation, the sum of all 
terms in the three-dimensional formalism must produce covariant results. In nuclear physics 
it is known that contributions which are of higher-order in the coupling are of increas- 
ingly shorter range. Such short-range contributions to hadronic interactions must always be 
treated in an essentially phenomenological manner. Therefore, it is expected that a trun- 
cation of the interaction at some finite order in the coupling will be useful for applications. 
It should always be appropriate to absorb the effects of the neglected higher-order graphs 
into phenomenological parameters, thus ameliorating the non-covariance of the theory. In 
particular, we show in Section [TTT] that the leading-order boost corrections to the interaction 



obtained by Forest, Pandharipande and Friar |27j are indeed contained within our truncated 
interaction. 

In Section we show that the crossed-box graph may be approximately rewritten as 
an iterate of the ladder kernel: 

KgaKWGcK®, (1.14) 

where the form of Gq is derived in Section [TV|. In particular, this approximation is exact in 
the high-energy (eikonal) and one-body limits, and may be systematically improved upon 
in other regimes. This leads us to propose the four- dimensional integral equation: 

T = K^(G + G C )T. (1.15) 

This equation has the correct one-body limit and results from an approximate resummation 
of the BSE kernel K. It reduces to the Wallace- Mandelzweig equation ( |1.11| ) if the depen- 
dence of on the time component of the relative four-momentum in the cm frame is 
neglected. A three-dimensional reduction of Eq. ( |1.15| ) using the method developed in Sec- 
tion H provides a systematic way to include retardation effects, and so improve on Eq. ( |1 . 1 1| ) . 
This reduction to a three-dimensional equation is performed in Section The dynamical 
boost of the three-dimensional interaction present in the equation is realized through simple 
dependence of the interaction on the total three momentum and energy of the system, thus 
eliminating the boost problem of Eq. ( |1 . 1 1| ) . 



In Section [VI] the predictions of the three-dimensional integral equation obtained in 
Section [V] are compared to those of five other bound-state equations: the ladder BSE 
and Eq. ( |1.15|) , both of which are, of course, four- dimensional equations, and the three- 
dimensional equations of Salpeter, Klein and Wallace & Mandelzweig. 

The results of Sections are presented concurrently for a scalar and a spinor field 

theory. In particular, in the scalar case we use the <j) 2 a field theory, defined by the Lagrangian 

£ = 7^0i<9 M 0i - m\4>\ + <9 M 02<9^2 - ml4>l + d^ad^a - \i 2 a 2 ) - g\m\$? - Qimtfjfco. 

(1.16) 
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Coupling terms include mass factors so that the couplings g\ and g% are dimensionless, thus 
ensuring that the limits m\ — > oo and m 2 — > oo correspond to the appropriate one-body 
limits. It is this field theory in which the numerical calculations of Section [VT| are performed. 
In the spinor case the Lagrangian is 

£ = ^(ij^ - mi)^i + ^2(27^ - m 2 )^2 - S^if^i - 92^20^2- 

The arguments of Sections PTT]^V| are, in fact, quite general, and, with small modifications, 
also apply to field theories involving other types of particles. 



II. A SYSTEMATIC PROCEDURE FOR THE REMOVAL OF THE 
RELATIVE-ENERGY DEGREE OF FREEDOM FROM A BOUND-STATE 

EQUATION 

Consider the Bethe-Salpeter equation for the bound-state vertex function T, i.e., 

r(p;P) = / 0^K(p,p';P)G o (p';P)T(p';P), (2.1) 

Here P is the total two-body four-momentum, and p and p' are the relative four-momenta. 
This equation is completely general, and applies to the bound-state vertex function in any 
field theory. Throughout this paper we find it convenient to omit the explicit integration 
from such equations, abbreviating them as follows 

T = KGqY. (2.2) 

Suppose that the driving term of this equation, K, is separated into two pieces: 

K = K l + K 2 , (2.3) 

where K\ is three-dimensional in the sense that it does not depend on the zeroth component 
of the relative four-momenta p and p' . For instance, if ( |2.1| ) was being solved in the two-body 
center-of-mass frame and K was a one-boson-exchange kernel then K% could be chosen to be 
the usual static one-boson-exchange kernel. Note that even though the whole driving term 
K is Lorentz covariant, because the restriction on K\ is frame dependent, the pieces K\ and 
K2 are not. The splitting of K leads to the following coupled equations for V: 

T = T l + T 2 , (2.4) 
T 1 = K 1 G (T 1 + T 2 ), (2.5) 

r 2 = i<: 2 G'o(r 1 + r 2 ). (2.6) 

Both these equations are four-dimensional, however, Ti has no dependence on the zeroth 
component of the relative four-momentum because of the defining condition of K\. Mean- 
while, the second equation may be formally solved to obtain: 

r 2 = [i-^ 2 G ]" 1 ^ 2 G'or 1 
= -r 1 + [i-^ 2 G ]" 1 r 1 . (2.7) 
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This last result may then be substituted into Eq. ( |2.5| ) to yield 



r\ = K x gr x , (2.8) 

where 

g = G + G K 2 g. (2.9) 

Eq. ( |2.8| ) then becomes a three-dimensional equation because the implied integrations over 
time-components of momenta only affect Q, i.e., it reduces to 

r 1 = K 1 (g)r 1 , (2.10) 

where 

(0> = J%$ev,,nn (2.1D 

Here we have used the same implied integration notation for the three-dimensional integral 
equation ( [2.10|) as for a four-dimensional equation. This practice continues below and the 
context should make it clear whether the equation in question is a three or four- dimensional 
one. 

Eqs. (|2.10|) and Q2.9|) are exactly equivalent to the original BSE ( |2.1|) . In order to 
reconstruct the Bethe-Salpeter amplitude T from the non-covariant three-dimensional Ti 
one must use Eq.( j2.7P rearranged into the form 

T = [1-K 2 G }- 1 T 1 . (2.12) 

The Green's function, (Q), defined by Eq. (|2.11|) is, in fact, the Fourier transform of the 
corresponding equal-time co-ordinate space Green's function Q{t', x' 1; x^; t, x 1; x 2 ). This may 
be seen by defining 

x^x^x^xa) = <5(a;°' - x 2 )5(xl - ^M^'n x' 2 ; x u x 2 ), (2.13) 

where g is the appropriate four- dimensional two-body Green's function, then inserting the 
integral representation of the delta functions and taking the Fourier transform . 

If K\ is now chosen to be a cm-frame instant interaction, K ms ^, then ( |2.9|) may be 
taken to define Q order by order in if - fT ms ^. This provides a systematic way to calculate 
retardation corrections to the Salpeter equation 



r 1 = K inst (G )r 1 , (2.14) 

which results from taking Q at zeroth order in K — K ms ^. In particular, at first order in 
K - A^ inst we find that: 

(S«> = (Go) + (G KG ) - (G )K inst (G ). (2.15) 

The use of this three-dimensional propagator in place of (Go) i n (|2.10| ) would therefore 
incorporate leading-order retardation corrections in the three-dimensional equation. 
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There is, however, a general way to account for all the relative-energy integrations in the 
modified Green's function (Q). Recall that the choice of Ki is subject only to the constraint 
that it should not depend on the zeroth component of the relative four-momentum in, say, 
the cm frame. Consider Eq. (|2.9|) rewritten as 

(Q) = (G ) + (G KG) - {Gq)Kx{Q), (2.16) 

and choose K\ such that (Q) = (G Q ). This provides the defining condition, 

K 1 = (G )- 1 (G Kg)(G )- 1 . (2.17) 

Since this choice of Ki means that (Q) = (G ) the full dynamics reduces to a three- 
dimensional integral equation with the free propagator (Go), i-e., 

r^i^Go)^, (2.18) 

and all the complexities of the relative-energy integrations are transferred to the interaction, 

It follows that K\ is the two-particle irreducible (2PI) interaction, where two-particle 
irreducibility is defined with respect to the three-dimensional propagator (G ). To show 
this formally, rearrange (GqKQ) as follows, 

(GoKQ) = (G KG (1 - KGo)- 1 ) + (G KG {[1-{K- K^G,}- 1 - [1 - ^Go]" 1 }) 

= (G TG ) - (G TG }K 1 (Q), (2.19) 

where T = K(l — GqK)^ 1 is the Bethe-Salpeter t-matrix defined by Eq. (|LT|). Use of 
Eq. ( tng| ) in Eq. ( PT7|) produces 

T\ = K x + Ti(G Q )Ki, (2.20) 

where T\ = (Gq)~ 1 (GoTGo)(Go)~ 1 is the three-dimensional t-matrix corresponding to the 
Bethe-Salpeter t-matrix T. Thus K\ is, indeed, the irreducible interaction which produces 
the t-matrix T\ using the propagator (Go). We therefore expect that the rules for the 
construction of K\ will be akin to those for the two-particle irreducible interaction in time- 
ordered perturbation theory. 



III. CONNECTION TO TIME-ORDERED PERTURBATION THEORY AND 

THE WORK OF KLEIN 

The connection to time-ordered perturbation theory (TOPT) emerges for a simplified 
dynamics in which only positive-energy states are kept in Go- For the case of a field theory 
of nucleons and pions it has been shown that the Green's function (Go + GqTGq) obeys the 
rules of TOPT in the no-anti-nucleon case pSySIfl . The same derivations suffice to prove 
the result for the scalar field theories defined in the Introduction. It therefore immediately 
follows that K\ is the 2PI TOPT amplitude in the case that only positive-energy states are 
kept in Gq. 
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We now provide a specific example of the derivation of the time-ordered perturbation 



theory amplitude by this route. Firstly, it is true in general that, if in Eq. ( 2.17 ) Q is replaced 
by G , then the three-dimensional kernel becomes 

K 1 = (G )- 1 (G KG )(G r l . (3.1) 

We note that this formula was given by Klein in his work on deriving three-dimensional scat- 
tering equations from four- dimensional ones [p3|-|2"5|j . Equations ( |2.17| ) and ( |2.9| ) may thus be 



thought of as providing a generalization of the formalism of Klein. Secondly, suppose that 
K is expanded according to the power of the coupling constant in each of its contributions, 



K = f^K {2i \ (3.2) 

i=l 

and for the moment only the second-order, or ladder, contribution to K is kept, i.e., ( p.l|) 
is rewritten as: 

K? = (G )- 1 (G K i2] Go)(G )- 1 . (3.3) 

We now calculate this amplitude in the absence of negative-energy states. First, the propa- 
gators are split according to Eq. (p..3[), and only the positive-energy pieces retained. Second, 
the inverse Fourier representations of all quantities are inserted, so that the time-ordering 
can be clearly elucidated. Third, a change of variables to time differences is made. Fourth, 
the relative-energy integrations are performed, generating delta functions on some of the 
time differences. Last, the integrals over these time differences are performed. This proce- 

(2) 

dure allows the calculation of K\ in any frame. Indeed, it is completely general, and can 
be used to calculate Ki no matter what K is chosen. 

If the total four-momentum of the two-body system in the frame of interest is written 
P = (E, P), then it is seen that K[ takes the standard form of TOPT: 

K? (p<, p; £1, P) = f*^ (- L_ + L_ ) , (3 .4) 

2a; V E+ — €\ — e 2 — w E + — e 2 — e\ — uj I 



where E + = E + ir], e\ = ei(pi), e 2 = e 2 (p2), u = + (p - p') 2 > and e[ (e' 2 ) has the 
argument p[ (p' 2 ) in place of pi (p 2 )- Here, mi and m 2 are the the masses of the two 
interacting particles, /i is the mass of the exchanged meson and the factor M. is defined by: 

_ f 4m 1 m 2 , for spin-zero particles, ,^ ^ 

1 1, for spin-half particles. 



Note that to extract the mass of the bound-state, M, one must set E = VM 2 + P 2 . The 
total and relative three-momenta are related to the individual particle momenta in the usual 
way: 

Pi + p 2 = Pi + P 2 = P; (3.6) 
zajPi - ^ip 2 = p; z/ 2 p'i - ^ip' 2 = p'; (3.7) 



S 



where v% and v% are any two real numbers which obey 1/1 + 1/2 = 1- Figure [l| shows the 
time-ordered diagrams corresponding to Eq. (3.4). 

The bound-state equation corresponding to this interaction is then: 

(Go)-V; e, pMp' ; e, p) = / -0-k[ 2 \ p ', p; p)V>(p; ^ p), (3-8) 

where we have omitted the factor of (2ti) 3 5^(p' — p) from (Go), which integrates out of the 
equation trivially. Here: 

< G °><" g - p >- E+ ^-^ (3 ' 9) 

where Af = A^pi) and = A 2 f (p2). Now suppose that we have a set of free single- 
particle positive-energy momentum-space wave-functions, u, for both the spin-half and scalar 
theories. In the spin-half case these will be the spinors, and their normalization may be 
chosen such that: 

A+(p)= M (p)m( P ) (3.10) 

(where spin indices have been suppressed). In the scalar case Eq. ( |3.10|) may also be enforced, 
since the u's may be chosen to be merely numerical factors, i.e., we may define 

«(p) = u(p) = -=L=. (3.11) 



Expanding: 

V(p; E, P) = u 1 ( Pl ) M2 (p 2 )0(p; E, P), (3.12) 



(with an implicit summation on spin indices in the spinor case) and manipulating Eq. ( |3.8|) 
leads to: 

{E-e' 1 -e' 2 )<t>{p'- 1 E,T>) = J -^-R?\p' , p; E, P)0(p; E, P), (3.13) 

with: 

K?\p', p; E, P) = ^(p'J^p^fV, p; E, P)n 1 (p 1 )n 2 (p 2 ). (3.14) 



Note that the TOPT interaction which appears in the bound-state equation ( |3.13| ) 
changes with total three-momentum of the system. This provides the dynamical boost 
of the interaction and it is straightforward to show that it incorporates the perturbative 
boost correction discussed in Ref. [p?H . That is, in the equal-mass case, to second order in 
P we have (neglecting p 2 /m and p /2 /m terms): 

if{ 2) (q; P) = k[ 2 \ q] 0) - ^^{ 2) (q; 0) - -Lp ■ qP ■ V q /?f } (q; 0), (3.15) 
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where we have written the interaction as a function of the momentum of the exchanged 
j 



meson, q 



P- 



If it is true, as argued in Ref . |27j , that this leading-order boost correction is sufficiently 
accurate for studies of three-nucleon systems, then the boost effects incorporated in the 
TOPT result ( |3.4| ) should be more than adequate. Indeed, we believe such a boost to be 
preferable to that of Eq. (|3.15 ), since it involves operators which are bounded as P —>■ oo, 
whereas, at large values of P the perturbative boost correction ( |3.15| ) diverges. 

The identification of the kernel K\ with the usual time-ordered perturbation theory 
interaction is only correct in the absence of negative-energy states. With the full four- 
dimensional propagator, the form of obtained from Eq. ( |3.3| ) in the manner described 
above Eq. (|3^) is: 



(G ) (p' ; E, P)K[ 2) (p', p; E, P) (G ) (p; E, P 

gmM 




(3.16) 



where Af = A^(pj) and Af = Af(p' i ). Note that here, and throughout the rest of this 
paper, the irj prescriptions can be generated by understanding the masses mi, m 2 , and fi to 
have a small negative imaginary part. 

Now the three-dimensional propagator is 



<G )(p;£,P) 



A+A+ 



+ 



AfA; 



E + ir] - e\ - e 2 —E + ir] — e\ 



t-2 



(3.17) 



(where once again we have omitted the momentum-space delta-function for notational sim- 
plicity), which is the Salpeter propagator [30|. If the negative-energy piece is removed, this 



becomes the Blankenbecler-Sugar |L2| or Logunov-Tavkhelidze |29| propagator. Note that 
in the spin-half case this propagator does not have a unique inverse, consequently in that 
case the derivation of the bound-state equation analogous to Eq. ( |3.13|) usually assumes that 
the H — and — h pieces of (Go) -1 are zero. 



The graphs corresponding to the terms in Eq. ( 3.16 ) are shown in Figure^. Graphs (a)- 
(d) are generated by the first line of Eq. ( ft.lq 1 ), graphs (e)-(h) by the second line, and graphs 
(i) and (j) by the third line. The fourth to sixth lines of ( |3.16|) generate the same ten graphs, 
but with particles and anti-particles interchanged. Note that the graphs drawn in Figure pi 
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(2) 

do not represent actual physical processes, but rather contributions to the "potential" K\ . 
As observed by Klein, if negative-energy particles are present, the rules for the construction 



of this interaction differ from those of standard TOPT p3p4 |. 

In general, it can be shown that the rules for the construction of the full two-body equal- 
time Green's function, (G + G TG ), for the case of positive-energy initial and final-state 
particles, are, in either of the two field theories given in the Introduction: 

1. Draw all topologically distinct, two-to-two graphs, which do not contain vacuum- 
vacuum subdiagrams, remembering that different time-orderings contribute to different 
graphs. 

2. Construct the expression for each individual graph exactly as in TOPT, save that: 

(a) All one-particle energies are to be relativistic. 

(b) If in some intermediate state of the graph both the initial-state particles are 
present with exactly the same momenta as in their initial state, i.e., their state 
is unaltered from the beginning of the graph, then in the TOPT denominator 
corresponding to that state the quantity e± + ei must be replaced by E. 

(c) Similarly, if both final-state particles are present with exactly the same momenta 
as in their final state, i.e., they have undergone their last interaction but are still 
present in some intermediate state, then in the TOPT denominator corresponding 
to that state the quantity e[ + e' 2 must be replaced by E. 

If both these last two conditions are satisfied then the denominator for the relevant state 
is the TOPT denominator for that state, but with e[ + e' 2 + e\ + e 2 replaced by 2E. Note 
that these are not the rules for the amplitude obtained from the Ladder BSE. The Green's 
function defined by these rules is the ++ — > ++ piece of the full two-body equal-time 
Green's function of the field theory. 

The substitution of these Es in the intermediate-state denominators can be shown to 
be due to a crucial difference between the way initial and final-states are treated in the 
calculation of (GqTGq) and in TOPT. In TOPT the initial and final times in the diagram 
are taken to minus and plus infinity respectively, thus guaranteeing that no interaction 
takes place before the beginning (or after the end) of the propagation of the two-particle 
state. However, when {GqTGq) is calculated in our work the initial and final times are kept 
finite. Consequently the first (last) event in the diagram need not be the creation of the 
0102 or NiN 2 pair, it may be the creation out of (destruction into) the vacuum of a (fxpa 
or NNa state. (See for instance the second graph in Figure |2|.) The possibility of such 
an event happening before or after the 0102 or NiN 2 propagation in some contributions 
to (GqTGq) leads to denominators which differ from the TOPT ones in those pieces of 
(GqTGq). It follows that, if the pieces of the Hamiltonian of the theory which couple 
(jxfta or NNa states directly to the vacuum are removed from the Hamiltonian, then the 
modification of the TOPT denominators is not necessary, and the rules for calculation of 
the two-particle to two-particle part of the Green's function (G TG ) become exactly those 
for the calculation of the TOPT Green's function. This equivalence of (GqTGq) and the 
time-ordered perturbation theory Green's function in the absence of such terms from the 
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Hamiltonian was demonstrated for a field theory of nucleons and pions by Kvinikhidze and 



Blankleider 28 



As discussed in Section [TI|, the use of a static interaction in the Bethe-Salpeter equation 
leads to the Salpeter equation (|2.14j) . This equation is often used for the interactions of two 
relativistic particles. It is straightforward to show that the static limit of the potential given 
by Eq. (|3.16|) is the usual instant interaction in the positive-energy sector. On the other 

hand, Eq. ( |3.16| ) predicts that the coupling between the ++ and states is suppressed by 

retardation effects which cause it to be a factor uj/m (in the equal-mass case) smaller than 
the rest of the static interaction. This factor is zero in the static limit. Thus, for a causal 
interaction, the correct static limit of the equation ( |2.18j ) with the interaction ( |3.16| ) is the 



Breit equation, i.e., Eq. ( |2.14|) with the states omitted from the calculation. 



In order to recover the results of Klein for two-meson exchange interactions, and so make 
connection with the recent work on two-pion exchange of Rijken and Stoks |3^j33H and Lahiff 
and Afnan |33J, we evaluate Q to first order in K — K\ in Eq. ( |2.17| ). This produces 

K x = (G )-\(G KG ) + (G KG (K - K 1 )G ))(G )- 1 . (3.18) 



Now expanding K as in Eq. ( |3.2p and dropping terms of higher than fourth order in the 
coupling constant leads to 

K{ 4) = (G )- 1 f(Go^ (4) G ) + (GoK^GoK^Go))^)- 1 - K[ 2) (G )K?\ (3.19) 



When applied to particle-particle scattering Eq. (|3.19| ) gives exactly the same results 



obtained by Klein, except that where Klein removed the iterated second-order three- 
dimensional interaction by hand, here the formalism provides the subtraction naturally. 

The procedure developed in Section || thus provides the generalization of Klein's method. 
The interaction and vertex function defined by Eqs. ( [2. 17[ ), ( |2.9| ), and ( [2.18| ) have the full 
relative-energy dynamics of the Bethe-Salpeter equation (|2.1|) included in them. 



IV. A SIMPLE FOUR-DIMENSIONAL EQUATION WITH THE CORRECT 

ONE-BODY LIMIT 



Now let us return to the four- dimensional BSE ( |Llp . The simplest BSE with the correct 
one-body limit is 



KxG T, 



(4.1) 



where Kx is the sum of all 2PI ladder and crossed-ladder two-to-two graphs The 
bound-state masses predicted by such an equation have recently been obtained by Monte 
Carlo integration in the Feynman-Schwinger representation of both scalar 3 field theory and 



scalar QED by Nieuwenhuis, Tjon and Simonov [35—39 ] . In general though, the kernel Kx 



is too complicated for ( |4.1| ) to be solved by standard means. Therefore, in this section we 
shall derive a four- dimensional equation which has a simple kernel, the appropriate meson- 
production thresholds and the correct one-body and high-energy (or eikonal) limits. A form 
of the following derivation appeared in HIT 
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Although it is impossible to rewrite a crossed graph exactly as an iterate of the ladder 
kernel , the leading contributions of these graphs to the high-energy and one-body limits 
are iterative. This can be easily shown in the case of the fourth-order crossed-box graph. In 
both of the field theories defined in the Introduction, this graph corresponds to an expression: 

Kf{k' 1 ,k' 2 -k l ,h) = 

/cftv 1 1 

(2^ (g - P2? - ~ + H ~ Vi) (pz _ h? _ ■ (4-2) 

where the lines have been assigned the momenta shown in Figure |3], and P is the total 
momentum, which is conserved: 

P = k 1 + k 2 = k[ + k' 2 . (4.3) 



The propagators di were given in Eq. (|1.3|) , and the factor M. was defined in Eq. ( p.5|) . 

Following Wallace and Mandelzweig fIS|]il]fl we define, for any four-vector q, quantities 
q\\ and q± , as follows: 

% = q\\h + Q±^ (4-4) 

where P is the unit four-vector in the direction of P and 

q\\=q- P. (4.5) 

Consequently in Eq. (|4.2|) the argument of the function d,2 may be rewritten: 

(k 2]l + k' n - p 2{] )P + k 21 _ + k' 2± - p 2± . (4.6) 

Now suppose that the momentum of particle 2 is large. This may occur because m 2 ^> rri\ 
(one-body limit) or because particle 2 has very high energy (eikonal limit). In either case its 
intermediate and final-state momenta will be largely unaffected by the presence of particle 
1, and so we may approximate the perpendicular components as unchanging, 

k 2± + k' 2± w 2p 2± . (4.7) 

Indeed, making the replacement (|4.7|) in ( f4.2| ) will not affect the value of K^ 1 in the limit 
m 2 — > oo. 

This argument shows that K% may be approximately rewritten as: 
K^(k' 2 ,k 2 ;P)^ 

1 1 7^^ (2) (^-P2)rf 1 (P-p 2 )4((A: 2 ||+^ || -p 2 ||)P + p 2± )^ 2 )(p 2 -A ;2 ), (4.8) 

where is the ladder BSE driving term: 

K*\q) = M 2 9l9 2 \. ■ (4.9) 
q z — fi z + ir] 

In operator notation: 
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J\ x ~ 



(4.10) 



with 

Gc^P! £2) = *(2tt)V 4 V - p)d x {p)d 2 {{E 2 -E,+ 
Here £7i and £3 are defined via: 

2£ 2 = fey + k' n ; 2E 1 = k n +k' n - 



)P-P±) 



(4.11) 



(4.12) 



and so in the cm frame E = E\ + E 2 . The propagator Gc defined by ([4.1 1|) therefore 
depends on the parallel components of the external momenta. Thus the use of operator 
notation in Eq. ( 4.10|) is not strictly correct. In Eq. ( 4.27 ) we redefine Gc in order to remove 
this dependence on external momenta. However, that change in Gc modifies the analytic 
structure of the amplitude defined by the corresponding integral equation. Therefore, for 
the present we persist with Gc defined by ( (4.11|) , and use an improper operator notation. 
Note that if the particles are on-shell in their initial and final states then: 



E\ = E l 



on 



E 2 + m\ 



rrin 



2E 



E, 



E™ 



E 2 + m 2 ,— m\ 
2E 



(4.13) 



These arguments show that Eq. (|4.10 ) will be exact in the infinite-mass and high-energy 
limits, thus demonstrating that the pieces of K% which survive in these two limits may, 
indeed, be written as iterates of K^ 2 \ 
Now suppose that Kx is written as 



K x = V + VG G K X . 



(4.14) 



In principle this is always possible, as ( |4.14| ) may be taken as a definition of V. At second 
order in the coupling we clearly have: 



vV> = kV> = kQ\ 

while the above argument shows that with this 

K { $ = V^G C V (2 \ 



(4.15) 



(4.16) 

: KM. 



in the high-energy and infinite-mass limits. Thus, a reasonable choice for V is V 
Equation (|4.14 ) then defines corrections to this choice. 

Once this V is chosen, Eq. ( |4. 14 ) and the BSE Q4.1| ) may be combined to yield an 
"improved" ladder BSE, which, in the two-body cm frame, after a change of variables to 
total and relative four-momenta, takes the form 



(2k) 



d 2 (E° n - Po , -p) + d 2 (E^ -p' + Po , -p) r(po, p; s). 



-ion 



(4.17) 



We stress that what has been done here is to take certain pieces of the Bethe-Salpeter 
kernel K x and rewrite them in the form K^GcK^ 2 \ K^GcK^GcK^ 2 \ etc. Consequently, 
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Eq. ( [4.171 ) is equivalent to a Bethe-Salpeter equation in which graphs other than one-meson 
exchange are approximately included in the kernel. Thus we expect that the solution of this 
equation may provide a better description of the dynamics of two-particle systems than the 
ladder BSE amplitude. 

This will be especially true for systems with one particle much heavier than the other. 
In these systems the one-body limit constitutes an important piece of the dynamics. Unlike 
the ladder BSE, Eq. (|4.17[ ) has the correct one-body limit. This may be shown as follows. 
For the on-shell vertex function p' Q = 0. If p' = the sum of the two particle 2 propagators, 
multiplied by M. is, in the spin-zero case: 



M 



d 2 (E» lL - p , -p) + d 2 (E 2 JlL + po, -p 



Am\m 2 



(E™ -p,y-el + i V (E™ + V of - e|. + i V 1 

while for the spin-half case the same combination becomes: 

(E 2 on - p ) 7 ° + 7 2 ■ p + m 2 (£? n +p )7° + 72-p + m 2 
(£on_ po)2 _ e 2 + i7? + (E§*+po)*-4 + iTi ' 

In the m 2 — > oo limit, E^ 11 — > e 2 — >• m 2 , and so this expression reduces to 

' ' 2mi = -2vri5(p u )2mi, 



Po + iri Po- %T) 
for the spin-zero case and 

- 27T^(p )A 2 + o 

for the spin-half case, where 

1 + 7o 



A 



2oc 



(4.18) 



(4.19) 



(4.20) 



is the positive-energy projection operator for an infinitely massive particle two. Thus, in the 
infinite m 2 limit Eq. (|4.17 ) yields the Klein-Gordon or Dirac equation for the wave function 
of particle 1 moving in the static a field generated by particle 2 



E 2 



J2 



or 



2m! 



(£ 7 ° - p' ■ 7 - mOVCp') 



d 3 p 



9i92 



(27r)s (i/ - p) 2 - a* 2 



^(p), 



d 3 p 



gm 



(2vr)3 ( P ' - P )2 



rV>(p), 



(4.21) 



(4.22) 



where 



^(P) 



E 2 — m 2 — p 



rr(0,p;£) 



(4.23) 



or 
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v>(p) 



1 



E'j — p • 7 — mi 



(4.24) 



and £ = E° n . 

Equation (|4.17|) also has the appropriate meson-production thresholds, and therefore 
could be used as a basis for four-dimensional calculations of two-body bound-state properties. 
However, the equation cannot be written in the form 



r = K (2) GT, 



(4.25) 



and consequently the method of Section [TT] cannot be applied to it. Therefore, we now seek 
an approximate version of Eq. ( [4.1 7|) which can be written in the form ( [4.25|) . 

In the on-shell vertex function p' Q — 0. In what follows we use p' — in the integrand 
also when the amplitude is not on-shell. This approximation provides a four- dimensional 
version of the three-dimensional Wallace-Mandelzweig equation obtained in Ref. [pHT9|j . In 
the center-of-mass frame it is 



IV; 0°*, E™) = i J Ak®(j/ - p)d 1 (Er + Po, P) 



(27T) 



d 2 (E? n 



2 - P0: - p ) + d 2 (E° n +p ,-p) T(p;E° n ,E° n ) 



Redefining Gc to be, in the cm frame, 

G c (p',p;E° n ,E™) = i (2n)^(p'-p)d 1 (E™+p ,p)d 2 (E° n +p , 



allows Eq. ( 4.26 ) to be written: 

r = K ( - 2 \G + G c )r. 

The following points are worth noting: 



(4.26) 



(4.27) 



(4.28) 



1. Because the new choice for Gc does not affect the on-shell fourth-order piece of the 
amplitude Eq. ( |4.10| ) is still exactly true on-shell in the one-body and eikonal limits. 

2. In the limit m 2 — > oo Eq. ( |4.26|) has the correct one-body limit, as can be seen by a 
similar argument to that given for Eq. (|4.17|) . 

3. If a pinch analysis of the singularities of the amplitude T defined by the scattering 
equation corresponding to ( |4.26| ) is performed then it is found that at order g 4 there is 
a production cut in the p' Q plane of the half-off- shell amplitude which extends upwards 
from 



Po = m 2 + fi- E% n , 



(4.29) 



and lies infinitessimally below the real axis. In the on-shell amplitude at sixth order 
this cut appears in the integrand and overlaps with the pole of the Green's function 
at p'q = -E° n — e 2 (q), so producing a cut in the _E° n plane of T which extends upwards 
from 
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= m 2 + ^. (4.30) 

This is not a particle-production threshold which exists in the full Bethe-Salpeter 
scattering amplitude. Hence its existence in the amplitude T defined by the scattering 
equation corresponding to Eq. (|4.26|) must be regarded as a deficiency of that equation. 



It arises because the method used to derive Eq. (|4.26Q is only guaranteed to produce 



the correct cut-structure for the second and fourth-order on-shell amplitudes. It should 
be noted that the vertex function given by Eq. ( |4.17| ) does not contain this anomalous 
threshold. 



The value of E at which the threshold ( |4.30| ) occurs in the on-shell amplitude is 



E = m 2 + | + Jm\ + m 2 fi + l r- In fact, since Equation ( |4.17|) was designed for use 



in the m 2 > m\ regime this is actually above E — m x + m 2 + fi- Therefore the theory 
has the correct threshold structure 

for all energies E such that m 2 — mi < E < mi + m 2 + \l. Furthermore, if mi = m 2 
then the threshold ( 4.30Q actually lies at E = mi + m 2 + /i, which is where the first 



production threshold of the crossed-box graph should be. Consequently, we think of 
this cut as representing the usual single-pion-production threshold of the crossed-box 
graph, but somewhat displaced if m 2 > mi. 

4. This equation was derived for particles of different mass. If we desire an equation for 
identical particles a symmetric form of the propagator Gc must be used. (We shall 
return to this point in the next section.) 

5. The above derivation could equally well be pursued for the exchange of vector particles. 
However, in, for instance, scalar QED, account must be taken of the seagull graphs. 
Furthermore, in the case of spinor QED more care must be taken, since an additional 
piece of the interaction may be generated when one attempts to write X « VGqV. 



(The interested reader may consult [19j for details on these points 



V. A THREE-DIMENSIONAL EQUATION WITH RETARDATIONS AND THE 

CORRECT ONE-BODY LIMIT 

We now apply the relative-energy integration method of Section [II] to the four- 
dimensional equation Q4.26| ). By so doing we first recover Wallace and Mandelzweig's original 
quasi-potential equation, and then calculate a three-dimensional kernel which includes first- 
order retardation corrections to this result. 

Observe that once Gc is defined by Eq. ( |4.27| ), Eqs. ( [4.14D and ( |4.1| ) are equivalent to: 



T = V(G + Gc)T, (5.1) 

where V is regarded as being defined by ( |4.14j ). This provides a method for calculating 
corrections to Eq. (|4.28| ), which was obtained by assuming that V = K {2 \ The result of 



applying the method of Section [TI] to Eq. ( |5.1[ ) is: 

Ti = V x {g)T x , (5.2) 
G = G + G c + (G + G c ) (V-Vi)G. (5.3) 
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If V is taken to be V^.j ng ^, where V^.j ns ^ is the static one-sigma exchange potential 

9i9%M 
; P -p') 2 + /i 



^inst(P >P) = ~ f „ ^2>„2 > ( 5 - 4 ) 



these equations reduce to 

Ti = V ainst G ET T u (5.5) 

which is a three-dimensional equation with a free two-body propagator that in the cm frame 
takes the form: 



G ET (p- E° n , £ 2 on ) = (Go + G c ) (5.6) 

^2 



A+A+ A+A7 



+ 



E? n + E§ n -ei-e 2 + iv E? n - E§ n -ei-e 2 + ir] 
A^A+ Af A 2 



' -E° n + E° n - ei -e 2 + i7 1 + -£° n - £ 2 on - e t - e 2 + irj ' (5 ' ?) 

where once again a factor of (27r) 3 <r 3 '(p' — p) has been omitted for notational simplicity. Note 
that this propagator is equivalent to that derived by Cooper and Jennings by a different 
technique Note also that in the case mi = m 2 this same propagator is referred to as 
the equal-time propagator by Tjon and collaborators |15 , 37| , 38| , ^I , ^2 1 . 

The desired generalization of Eq. ( |5.5| ) which incorporates retardation and boost effects 
follows from applying the ideas of Sections |I| and |T| to Eq. (|5.3| ). If V% is chosen such that 
the last term in Eq. ( [5.3| ) is zero, the formula 

V x = G^T ((G + G c )Vg)G ET , (5.8) 

is obtained for the general three-dimensional interaction to be used in 

Fx = VxG ET Y x . (5.9) 

A first step in the inclusion of retardation and boost effects may be taken by following 
Klein's work and replacing Q by Go + Gc- This yields what we refer to as the first-order 
relative-energy integration result: 

Vx = G ET AG ET (5.10) 
A=((G + Gc)V(G + Gc)). (5.11) 

We now take V equal to and calculate A in the center-of-mass frame. In order to 
do this we first decompose the cm frame propagators: 

*(p) = E — lAE tt; * = 1 > 2 > ( 5 - 12 ) 

where pi may take on the values ±1. It follows that the amplitude A may be written as 

A= £ h£'i$'A(p' x p 2 ^pxp 2 )i^I^. (5.13) 

Plp' 1 P2P' 2 
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The sixteen contributions A(p' 1 p' 2 <— p\p%) may be calculated by the same method used 
to obtain the Klein potential (|3.16|) . They are found to take on a slightly different form 
depending on whether p\p\ is equal to plus or minus one. In order to simplify the expressions 
for the A' 8 we define the following quantities: 



9192M 



2uj 



on 



r! 77 ,( 



on 



If p[pi = 1 then: 



A(p[p' 2 <- P1P2W, p; ^ on , E™) = F 

+ F 

Conversely, if p\p\ = — 1 then: 



e[ + e' 2 t\ + e 2 — uj 
1 1 



1,2. 



1 



(5.14) 
(5.15) 



+ 



1 



ei + e 2 e 2 + e' 2 - uj 
1 1 



ei + e' 2 e 2 + e' 2 - u 



ei + e 2 - u; ei + e 2 



(5.16) 



^(p / lP , 2-PlP2)(p',p;i? 1 on ,i? 2 on, 



= F 
+ F 



1 



1 



e% + — uj ei + e 2 — u; 
1 1 



1 



1 



+ 



ei + e' 2 e 2 + e' 2 - lu 



+ 



ei + e 2 
1 



e 2 + e 2 — uj 
1 



ei + e 2 - uj e 1 + ei - lu 



ei + e' 2 



ei + ei — e 2 + e 2 — uj 



ei + e 2 



(5.17) 



We note the following points about the interaction defined by Eqs. (|5.7|) , (|5.10| ), ( ^.13|) 
( CT ) and ( gT7p . 



1. As expected from the four- dimensional equation ( |4.26| ) it assumes a static form in the 
limit m 2 — > 00. In other words, it has the correct one-body limit. (Recall that we 
assumed that particle 2 was the heavier of the two.) 



2. In the static limit of the interaction defined by Eqs. (|5.16 ) and ( 5.17 ) all couplings 

to the states are suppressed. On the other hand, the other couplings tend to the 

same static limit as the +H > ++ piece of the interaction. Thus, in this limit the 

+H — > H — coupling, which is due to Z-graphs, is correctly given by the instantaneous 

exchange interaction, and only the states need be omitted from a calculation. The 

role of the H — and states in this limit is therefore consistent with the one-body 

limit which motivated Eq. (|5.9|). 



3. Also as expected from (|4.26| ) A(++ < h+) has a singularity structure which is differ- 
ent to that of the ordinary time-ordered perturbation theory amplitude: it contains a 
cut beginning at E 2 = m 2 + 1 . The same comments made about the unusual singularity 
structure of ( |4.26| ) apply to this cut. 

4. Because we have used a form of Gc derived by considering the limit m 2 — > 00, the 
interaction is not symmetric under the interchange of particle 1 and 2 labels. 
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Points (2) and (3) are not a reflection of the underlying physics of the meson-exchange 
process, but rather of the particular iterative form we used in our attempt to sum some of 
the higher-order graphs in the kernel of the BSE (O). 



The choice of Gc used above involves an approximation which is appropriate in the regime 
m 2 > mi. For equally massive particles a more appropriate choice for the propagator Gc is: 



G c (p>,p;E° n ,E° n ) = l {2^\p'-p)- dx(E™ + p , p)d 2 (E™ + p , 



d 1 (E^-p ,p)d 2 (E^-p ,-p 



(5.18) 



Note that this choice is label symmetric, and hence the resulting bound-state equation may 
be used for identical particles. The Gc defined by ( |4.27| ) did not have this property. 

When Eqs. ( |5.10|) and (|5.11|) are applied with this symmetrized form of Gc the result is 
a V\ of the form ( |5.13| ), but with an A(p[p' 2 <— P1P2) which may be written in any p-spin 
channel as: 



A {P,P,Ex ,E 2 



F 
~2 



+ 



1 



+ 



1 



e 'i + e 2 \ e i + e 2 ~ u e i + e 2 — uo ei + e[ — oj 
11 1 11 



+ 



e 2 + e' 2 - U / 
1 



e[ + e 2 e\ + e' 2 — to e 2 + e' 2 — uo 
1 1 1 



+ 



e[ + e 2 e[ + e 2 — to t\ + e[ — uo 
1 1 1 



ei + e[ - uo ei + e' 2 - uj e x + e 2 e 2 + e 2 - u e[ + e 2 - u e x + e 2 



+ 



1 



e 2 + e 2 



e x + e 2 



+ 



e 2 



ex + e\ 



OJ 



with all symbols defined as above. 



ei + e 2 



(5.19) 



VI. COMPARISON OF THE BOUND-STATE PREDICTIONS OF DIFFERENT 
THREE AND FOUR-DIMENSIONAL INTEGRAL EQUATIONS 



Having derived the interaction which is to be used in Eq. ( |5.9|) we may now compare the 
two-body bound-state properties predicted by this integral equation with those properties 
obtained from other three and four- dimensional calculations. This is done in the scalar field 
theory defined in the Introduction. The masses m 1; m 2 , and p are chosen to be mx = m 2 = m 
and p = 0.15m. Units are then chosen so that H = c = m = 1. 

The three-dimensional equations to be considered are: 

1. The Salpeter equation, which may be written: 

rs = V ainst (G )T s , (6.1) 



with V^^- given by Eq. (|5.4|) and (G ) by Eq. (|3.17|) . Note that in this scalar field 



theory Eq. ( |6.1| ) is exactly equivalent to the Blankenbecler-Sugar equation obtained 
from the ladder BSE. 
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2. The Klein equation for a one-sigma-exchange interaction, which is Eq. ( |2.18| ) with the 
interaction defined by Eqs. (|3.16|) and ( ft.17 ). 



3. The equation proposed by Wallace and Mandelzweig, which we refer to as the ET 
equation: 

Vet = V^j^GetTet, (6.2) 

where Get was given in Eq. (|5.7|). 

4. The retarded ET equation, Eq. (|5.2|) , with the label-symmetric interaction obtained 
in the previous section, i.e., with V\ defined by (|5.19|) , ( |5.13| ), ( |5.10| ) and 



The four-dimensional equations these four calculations are to be compared to are: 

1. The ladder Bethe-Salpeter equation. 

2. Equation (|4.26|) which, unlike the ladder BSE, has the correct one-body limit. 



We consider the six bound-state equations as falling into two groups: the first contains 
the ladder BSE and the two three-dimensional reductions thereof: the Salpeter equation and 
the Klein equation. The second group contains the four- dimensional Wallace-Mandelzweig 
equation ( |4.26| ) and the two three-dimensional equations derived from it: the equal-time 



equation ( |5.5|) and the first-order relative-energy integration equation (|5.9|) with the inter- 



action Vi defined by Eqs. ( |5.19| ), (|5.13|) , ( |5.10| ) and ( |5.7| ). Within these two groups the three 



equations may then be thought of as being zeroth (instant potential equations), first (Klein- 
type) and infinite (four-dimensional equations) order relative-energy integration results. 
At any given bound-state mass the equation: 

r(m) = M(m)T(m), (6.3) 

whether it be a three or four- dimensional integral equation, is an eigenvalue problem. The 
coupling constant gig 2 which appears in the kernel must be chosen such that M(m) has 
eigenvalue one. If we are searching for the ground-state of the system then we may assume 
that T is an S-waye state, thus implying that the function T has no angular dependence. The 
integration over p may then be easily performed. The resultant one (or two) dimensional 
integral equation may then be discretized, reducing the problem to one of solving for the 
eigenvalues of a matrix version of M. Note that in the four-dimensional equations it is 
necessary to perform a Wick rotation in the variable p$ before the kernel M is discretized, 
as otherwise the analytic structure in the kernel makes this procedure numerically unstable. 

When this approach is implemented for the equations listed above the six bound-state 
spectra shown in Figure [| are obtained. The following observations may immediately be 
made: 

1. The two instant formalisms (Salpeter and ET) predict more binding at a given cou- 
pling than either their corresponding four-dimensional equation or their corresponding 
Klein-type interaction. 
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2. The calculations in the second group always predict deeper binding than the corre- 
sponding calculations in the first group. In other words, the inclusion of pieces of 
higher-order graphs in the kernel of the integral equation yields deeper binding at a 
given coupling. 

3. The first-order relative-energy integration (Klein- type) interactions give predictions 
for the bound-state spectra which are significantly closer to those of the full four- 
dimensional calculation than the spectra of the instant equations. 

All of these results can be understood on physical grounds. Point 1 arises because the 
inclusion of retardation in the scalar field theory always reduces the amount of attraction 
in the bound-state. An example of this is provided by one-sigma exchange. When the 
full retardation corrections predicted by time-ordered perturbation theory are included the 
instant potential (|5.4j) is replaced by the retarded interaction ( |3.4j ) (provided the effects of 
negative-energy states are ignored). If E < mi + ni2 this will always be a negative number 
that is smaller in magnitude than the instant potential. Meanwhile, the inclusion of pieces 
of additional higher-order graphs leads to deeper binding since in the <p 2 a theory all time- 
ordered perturbation theory graphs are attractive in the bound-state region. Most significant 
though, in our opinion, is point 3, which shows that the first-order relative-energy integration 
(Klein-type) potentials do a better job than the instant interactions of approximating the 
four-dimensional equation which both are derived from. 

All six of these calculations aim to sum the main contributions to the two-to-two am- 
plitude. However, Figure |] shows that they give different results, and so they cannot all be 
correct. The question therefore arises: what is the binding energy at a given coupling if the 
full scattering series is summed? A partial answer to this question was recently provided 
by Nieuwenhuis and Tjon, who summed the series of all ladder and crossed-ladder diagrams 
using Monte Carlo integration techniques in the Feynman-Schwinger representation, and so 
obtained two-body bound-state masses in this scalar theory |35|-|38|j . Their results, calcu- 



lated at six different coupling strengths, show that the sum of all such graphs predicts deeper 
binding than any of our six integral equation calculations. In particular, the ET appears 
to do the best job of reproducing the bound-state spectrum obtained by Nieuwenhuis and 
Tjon. At first sight this appears to recommend the ET as the best formalism for doing 
three-dimensional calculations. However, two points must be remembered before a definite 
conclusion is drawn. 

Firstly, the reason that the Nieuwenhuis- Tjon calculation predicts so much more binding 
than the two four- dimensional integral equation calculations performed here is precisely the 
fact discussed in connection with point 2 above: in the 2 cr theory all time-ordered per- 
turbation theory graphs give attraction in the bound-state region. Hence the Nieuwenhuis- 
Tjon calculation, which includes many more such graphs than both the ladder BSE and 
four- dimensional Wallace-Mandelzweig equation calculations, must predict more binding. 
Therefore, we claim that the Salpeter and ET equations do a "better" job of imitating 
the full Nieuwenhuis- Tjon calculation than the ladder BSE and four-dimensional Wallace- 
Mandelzweig equation because instant calculations such as the Salpeter and ET ones ignore 
two competing effects. 

1. Higher-order graphs are left out of their kernel. If included these graphs would lead 
to deeper binding at a given coupling. 
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2. Exchanged-meson retardation is completely ignored. If included this retardation would 
lead to shallower binding at a given coupling. 

In other words, we claim that the ET formalism's apparent success is due to the cancellation 
of these two neglected effects and cannot necessarily be interpreted as a signal of that 
approach having the "right physics" for the problem. 

Secondly, none of these calculations have been performed in the way that calculations 
in few-body hadronic physics are usually performed. When such bound-state calculations 
are done in hadronic physics one usually acknowledges that the use of any integral equation 
necessarily means the neglect of some graphs of the theory, and so one regards the couplings 
that appear in that equation as effective. These couplings are then usually fixed by fitting 
to some set of observables, and then certain different observables are predicted. This makes 
the connection to an underlying field theory somewhat more tenuous, but it does allow 
calculation to proceed in situations where: 

• The value of the coupling which appears in the "fundamental" Lagrangian is not 
known. 

• The problem is highly non-perturbative and consequently any truncation is somewhat 
suspect. 

Therefore, to test the procedure of Section || for deriving three-dimensional integral 
equations in a way consistent with that in which it would be used in hadronic physics, we 
fix the mass of the bound-state at M — 1.95 and calculate the wave-functions 

*P(p';E)=g{p';E)r(p';E), (6.4) 

for each of the four different three-dimensional calculations. (Here g is the free two-particle 
Green's function for the relevant calculation.) Each three-dimensional wave-function must 
be normalized according to 

/ T^r^ (p ' ; jE) ^" (p '' p; E)ij{p] E) = 2E > (6 - 5) 

where G is the full two-particle Green's function appearing in the three-dimensional equation 
from which T was obtained. These wave-functions are then compared to the four- dimensional 
wave functions, integrated over the zeroth component of the relative four-momentum: 

W\E) = -Ly dp' g(p^E)T(p^E). (6.6) 

Here g is the appropriate free two-particle four- dimensional Green's function. The four- 
dimensional vertex function used in ( |6.6|) is normalized according to: 

/ ^tHp^^ ^ (6 ' 7) 

where G is the full two-particle Green's function appearing in the four-dimensional equation 
from which T was obtained. The left-hand sides of Eqs. ( |6.5|) and ( |6.7|) may be reduced 
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to two (and four)-dimensional integrals using the fact that the S'-wave wave-function if> is 
angle-independent. Note that in the four-dimensional case, Wick rotation in the zeroth 
component allows use of the vertex function on the po or p' Q imaginary axis. 

When the normalized wave-functions of the equations in group 1 (group 2) are com- 
pared the results shown in Figure |5] (|6p are obtained. It is immediately seen that all the 
wave-functions at this bound-state mass have essentially the same features. However, in 
both figures it is clear that the three-dimensional equation with the first-order relative- 
energy integration, or Klein-type, interaction reproduces the integrated four- dimensional 
wave function considerably more accurately than the instant calculation does. This shows 
that the procedure given in Section 2 does provide a way to systematically derive three- 
dimensional equations which are closer to their four-dimensional counterparts than those 
which are obtained from merely adopting an instant interaction. 



VII. CONCLUSION 

Quasi-potential equations (QPEs), despite their success as a basis for hadronic phe- 
nomenology, cannot be improved upon systematically. This is because whenever they are 
used beyond first-order they predict amplitudes with unphysical singularities. These sin- 
gularities also appear whenever attempts are made to boost the interaction appearing in 
the QPE from one frame to another. They arise from these equations' use of delta-function 
constraints on the relative four-momentum. 

Because of this shortcoming of QPEs, in Section || of this paper we sought, and found, 
a systematic procedure for deriving three-dimensional equations from four- dimensional ones 



which does not involve the use of delta functions. As shown in Section |T|, this procedure is 
akin to the work of Klein. It allows the derivation of a three-dimensional kernel which in- 
cludes, in a systematic expansion which may be pursued to any desired accuracy, the effects 
of the relative-energy integration present in the four- dimensional integral equation. If the 
procedure is applied to infinite order, a result equivalent to the original four- dimensional 
equation is obtained. At any order the kernel derived has a simple interpretation in terms 
of the diagrams of time-ordered perturbation theory, except that the presence of negative- 
energy states requires additional rules. The resultant three-dimensional equations therefore 
include more of the meson-retardation effects than quasi-potential equations, and incorpo- 
rate the effects of a dynamical boost. 

Since the BSE with any kernel which is a finite sum of Feynman graphs does not have 
the correct one-body limit, if the procedure of Section |H| is directly applied to any solvable 
BSE a three-dimensional equation without the correct one-body limit is found. However, 



as demonstrated by Wallace [[L8[], and recapitulated here in Section [TV|, a Bethe-Salpeter 
equation with the correct one-body limit may be reorganized so that the pieces of the kernel 
which contribute at leading order in the one-body limit take on an iterative form. Hence a 



four- dimensional equation, Eq. (|4.17p , with the correct one-body limit and meson-production 



thresholds is derived. A further approximation leads to a simple four-dimensional equation to 
which the procedure of Section |T|may be applied, Eq. (|4.26| ). In making that approximation 
some of the meson-production thresholds are displaced. Nevertheless, the equation still has 
exactly the correct cut-structure for mi — nil < E < mi + m 2 + [i. 
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In Section [V| the procedure of Section || was applied to the four-dimensional equation 
Q ). The result is a simple three-dimensional equation with the correct one-body limit 
which has straightforward boost properties. The main shortcoming of this equation is the 
fact that certain thresholds of the original Bethe-Salpeter wave function are modified in the 
three-dimensional wave function. However, this modification is due to the approximation 
in moving from (|4.17| ) to ( }4.26|) , rather than to any deficiencies of our procedure for deriv- 
ing three-dimensional equations. It could be systematically corrected for, at the price of 
complicating the three-dimensional kernel. 

Finally, in Section |VT| we compared and contrasted the bound-state properties predicted 
by six different bound-state equations in a scalar field theory. A first group of equations 
contained the ladder Bethe-Salpeter equation and two three-dimensional equations based 
on it: the Salpeter equation, and the equation obtained by applying our method to first 
order — which we referred to as the Klein equation, since the formula obtained for the three- 
dimensional interaction appears in Klein's work. In the second group were Eq. ([4.26|), the 



so-called equal-time equation and the equation derived in Section [VI|. It was found that 
in each group of three equations the equation derived by the Klein-like delta-function-free 
reduction technique of Section |H] approximated its parent four- dimensional equation better 
than the corresponding instant equation did. However, the bound-state spectra from the 
instant ET equation lies closest to the bound-state spectrum obtained from the sum of 
ladder and crossed-ladder graphs, as calculated by Nieuwenhuis and Tjon |35H38|. This can 



be understood on physical grounds and, we argued, is not necessarily a recommendation 
for the use of the ET equation in physical systems. Indeed, when the wave functions of a 
two-body system in a scalar field theory are examined it is seen that the equation derived 
in Section [V] does a much better job of reproducing the integrated four-dimensional wave 
function than the ET equation wave function does. 

The results of Section [VI| are in harmony with recent results for phase shifts in scalar- 
scalar scattering obtained by Lahiff and Afnan They found that the Klein method 
reproduces the phase shifts obtained from the corresponding Bethe-Salpeter equation much 
better than a Blankenbecler-Sugar calculation with an instant interaction. 

These formal and numerical results indicate that the procedure of Section [TT] is successful 
in its goal of providing a way to systematically obtain three-dimensional equations from 
four- dimensional ones without the use of delta functions. In particular, when applied to the 
four- dimensional equation (|4.26| ) this procedure yields a three-dimensional equation which, 
unlike the corresponding quasi-potential equation Q5.5|), has a well-defined boost. This would 
seem to make such an equation a good starting point for calculations in few-hadron systems. 
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FIGURES 



FIG. 1. The two graphs which contribute to the second-order three-dimensional kernel, if only 
positive-energy intermediate states are included in the calculation. 

FIG. 2. The ten graphs which contribute to the second-order three-dimensional kernel if propa- 
gation into negative-energy states is allowed. Note that these graphs can be interpreted as applying 
to positive or negative-energy particles. In fact, each graph shown contributes twice: the second 
time with positive-energy particles going forward in time interchanged with negative-energy parti- 
cles going backward. The dotted lines represent the initial and final times in each graph. Observe 
that in all but graphs (a) and (e) some interactions take place outside the interval [U,tf]. 

FIG. 3. The crossed-box graph, showing the momentum labels used in the text. 

FIG. 4. A plot of bound-state mass (^) versus coupling strength (^p) for the seven different 
calculations discussed in the text. The points with error bars are taken from and represent the 
FSR calculation of Nieuwenhuis and Tjon. The other calculations were all performed using three 
or four-dimensional integral equations, with the legend as indicated in the figure. 

FIG. 5. A plot of normalized wave-functions ktp{k) versus momentum k (in units such that 
m = 1), for a bound-state with M = 1.95m. The different wave functions shown are the integrated 
ladder BSE wave function (solid line); the result of the Klein potential calculation (long dashed 
line) and the Salpeter equation wave function (dotted line). 

FIG. 6. A plot of normalized wave-functions ki^{k) versus momentum k (in units such that 
m = 1), for a bound-state with M = 1.95m. The different wave functions shown are the 
integrated four-dimensional Wallace-Mandelzweig equation wave function (thick solid line); the 
wave-function resulting from the first-order relative-energy integration method applied to that 
equation (dash-dotted line) and the ET wave function (short dashed line). 
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